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We study the efficiency of the Atomic Frequency Comb storage protocol. We show that for a given optical 
depth, the preparation procedure can be optimize to significantly improve the retrieval. Our prediction is well 
supported by the experimental implementation of the protocol in a Tm 3+ :YAG crystal. We observe a net gain 
in efficiency from 10% to 17% by applying the optimized preparation procedure. In the perspective of high 
bandwidth storage, we investigate the protocol under different magnetic fields. We analyze the effect of the 
Zeeman and superhyperfine interaction. 
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I. INTRODUCTION 

The prospect of a quantum memory (QM) for light natu- 
rally emerges in the atomic physics community since it deals 
with the most fundamental interaction of light with matter. 
Atomic vapors have specifically drawn the attention because 
they are model systems with good coherence properties. Ex- 
perimental results rapidly follow with convincing proof-of- 
principle realizations. Most of them are based on stopped- 
light experiment and derived from the Electromagnetically In- 
duced Transparency phenomenon. More generally, the pro- 
tocols exploit a two-photon transition to benefit from a long 
live time in the ground state (spin coherence) [ 1]. A complete 
storage requires a large optical depth which also limits the 
signal bandwidth. As a consequence, the storage bandwidth 
is usually narrow (megahertz range in practice). The quest 
of wideband storage requires drastic changes. The protocols 
and the material should intrinsically offer a broad interaction 
linewidth. Inspired by impressive realizations in the classical 
domain J2, [H, the photon-echo techniques in rare-earth-ion 
doped crystals (REIC) appear as a prime alternative. 

REIC have indeed also remarkable coherence properties at 
low temperature. The large natural inhomogeneous broaden- 
ing gives a large potential storage bandwidth. It is not a limi- 
tation for the coherent transient phenomena. It has to be com- 
pared to the narrow homogeneous width (3-4 orders of mag- 
nitude narrower) to yield a high number of spectral channels 
10, Ht] . This is alternatively interpreted as a high multi-mode 
capacity [6], crucial for quantum repeaters application. 

Despite efficient demonstrations of classical storage 0], the 
photon-echo techniques should also be consider with precau- 
tion. A QM for light imposes different constraints. For ex- 
ample, the standard two-pulse photon echo (2PE) can be ex- 
tremely efficient benefiting from intrinsic amplification (pop- 
ulation inversion) of the medium. This regime is inappro- 
priate for a QM since it deteriorates the fidelity |8y. New 
specific protocols then emerge keeping the advantage of the 
technique (efficiency and bandwidth) but getting around the 
limitation of the 2PE (fidelity degradation): The Controlled 
Reversible Inhomogeneous Broadening (CRIB) |§, [Till 
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and Atomic Frequency Comb (AFC) lfl2ll . Theses techniques 
fully use the primary feature of the REIC (narrow homoge- 
neous width, large inhomogeneous broadening, spectral hole- 
burning preparation procedures ...). 

The CRIB protocol has been successfully implemented 
with record efficiencies |[l3l [Till . It does not suffer from the 
previously mentioned limitations of the 2PE because the pop- 
ulation mainly stays in the ground (no inversion). The CRIB 
has been implemented experimentally for weak light fields at 
the single photon level 01511 . This recent experiment shows 
a low noise level given by technical limitations. In the same 
context, the AFC has been recently proposed lfl2ll . Experi- 
mental results followed immediately afterwards demonstrat- 
ing also good noise feature and large efficiencies ITrjl [l7il . 
As compared to this realizations where the signal is stored 
directly in the optical coherence (optical AFC), the storage 
has also been performed in the spin coherences (Raman AFC) 
01811 . This critical step opens the way for long memory time, 
on-demand retrieval and high multi-mode capacity. 

The CRIB and the AFC have a lot in common since they 
are both inspired by the photon-echo technique. They both in- 
volve a rephasing of the optical coherences. The long storage 
time is obtained by a Raman transfer to the long-lived hyper- 
fine coherences. The main difference comes from the prepa- 
ration procedure. For the CRIB, a narrow absorption pro- 
file is tailored within the natural inhomogeneous broadening. 
An electric field gradient is then used to control and reverse 
the artificial broadening of a Stark effect sensitive medium. 
The hole-burning selection procedure has been optimized and 
may require a large number of independent pulses lfl9ll . For 
the AFC, the absorption profile exhibits a periodic distribu- 
tion of very absorbing peaks. The preparation procedure is 
based on spectral hole-burning towards a shelving state as in 
CRIB. A periodic optical pumping spectrum produces the cor- 
responding modulation of the population. This is sufficient 
to produce an optical AFC echo without external control of 
the atomic transition frequency. On the one hand, the ex- 
perimental implementation is not restricted to Stark sensitive 
REIC. On the other hand, a rudimentary preparation proce- 
dure based on the repetition of pulse pairs has allowed proof- 
of-principle demonstrations with relatively large efficiencies 
lfl6i[T7ll . This two points explain the rapid experimental pro- 
gresses lllq, I17I Il8n since the theoretical proposal of the AFC 
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In this paper, we study the preparation procedure. We 
specifically optimize the pulse sequence to improve the ef- 
ficiency of the optical AFC. We experimentally design the 
optical pumping spectrum to obtain the appropriate atomic 
comb. We also investigate the influence of the magnetic field. 
It changes the level structure and then the broadband optical 
pumping dynamics. This study is relevant in the large band- 
width regime required for quantum repeaters application. 

The paper is arranged as follows. In section [II] we show 
how the preparation sequence can be optimized for an efficient 
optical AFC retrieval. The results are analyzed by modelizing 
the population dynamics during the preparation stage. In sec- 
tion[Ill] we study the efficiency under a varying magnetic field 
and shows the influence of the Zeeman and superhyperfine ef- 
fects. 



II. COMB TAILORING AND STORAGE EFFICIENCY 

The generation of an echo out of a spectrally periodic 
medium is the key ingredient to explain the AFC storage. 
This is interpreted as a dipole rephasing in the time-domain 
1 12, Il7n or alternatively as diffraction on a spectral grating 
IU7l uOll . This latter approach is fruitful to calculate the re- 
trieval efficiency. We then describe the ideal comb shape that 
lead to the maximum efficiency. This analysis is well sup- 
ported the experimental results obtained in Tm 3+ :YAG . 



A. Optimized AFC shape 

The key idea of the AFC protocol is the storage capacity 
of a spectrally periodic medium. This observation connects 
the AFC to the well-know three-pulse photon echo (3PE) or 
more generally to space-and-time-domain holography. For the 
3PE, the first two engraving pulses create a population modu- 
lation. The third pulse is diffracted and produces an echo. The 
AFC is a direct descendant of the 3PE in that sense. To in- 
crease the population modulation depth by accumulation, the 
engraving sequence can be repeated. The AFC preparation 
proc edures used so far also involve accumulated pulse pairs 
lll6l[l7ll . The AFC has also some major differences with the 
3PE which makes it especially interesting for quantum stor- 
age. On one side, the protocol implies a Raman transfer to 
obtain a long memory time, on-demand retrieval and perfect 
efficiency. On the other side, the comb shape is a very sin- 
gular periodic modulation. In the limit of perfect retrieval, 
narrow absorbing peaks should be separated by fully transpar- 
ent regions [ 12]. This situation requires a large initial optical 
depth. This is usually a limiting factor in practice. We now 
can see how to optimize the comb shape to obtain the largest 
efficiency for a given optical thickness. 



the diffraction of a spectral grating. This model has been ex- 
plained previously, we here simply review the main results 
MM- 

We use the Bloch-Maxwell equations assuming the slowly 
varying amplitude and the rotating wave approximations I2H1 . 
In the weak signal limit ll22tl . the field propagation and atomic 
evolution are given by 



• f a(A)P(A;z,t) 



d z n{z,t) + \d t n{z,t) 

2ttJa (1) 
d t V(A; z,t) = - (iA + 7) V(A; z, t) - iCl(z, t) 

where t) is the Rabi frequency and V(A; z, t) = 
u(A;z,t) + iv(A;z,t) is the polarization including the in- 
phase and out-of-phase components of the Bloch vector (7 is 
the homogeneous line- width). We model the inhomogeneous 
broadening by a frequency dependent absorption coefficient 
a (A) including the coupling with individual dipoles and the 
inhomogeneous lineshape. 

The spectral periodic structure produces a series of echos. 
We decompose the absorption coefficient distribution in a 
Fourier series a (A) = a n e~ mAT (1/T is the comb pe- 
ri 

riod). Successive echos whose respective amplitudes are a p 
appear at times pT i.e. Vt(z, t) = a p (z) SI (0, t — pT), 

p>0 

p > for causality reason. 

We'd like to forewarn the reader. We here summarize the 
propagation through the spectrally prepared medium. For rea- 
son of simplicity, the comb is modeled by a frequency selec- 
tive absorption coefficient. It may give the wrong impression 
that the retrieval is a pure absorption effect. The causality im- 
poses that the successive echos appear after the arrival of the 
signal. This natural assumption is actually strong because it 
connects the field amplitude and phase by the Kramers-Kronig 
relation. The interplay between the field amplitude, phase and 
energy will be specifically addressed later in section Hi A 31 

So far, the different amplitudes of successive echos are 
recursively calculated from the propagation equation in the 
Fourier space. We here focus on the amplitude of the first 
echo. It defines the retrieval efficiency \ai (L) | 2 (L is the 
crystal length) and is given by successive integration with the 
boundary conditions ao (0) = 1 matching the incoming field 
intensity and a\ (0) = 0: 



d z a Q (z) 



1 



a a Q (z) 



(2) 



d z ai (z) = -- [a ai (z) + 2a_ia (z)] 



We finally obtain an analytic expression for the efficiency con- 
necting the comb shape and the first two coefficients of the 
absorption coefficient Fourier expansion. 



1. Retrieval efficiency 

We here calculate the retrieval efficiency. In the spectral 
domain, it is deduced from the comb structure. We consider 



7] (L) = \a-iL\ e 



-a L 



(3) 



It has been independently calculated by Sonajalg et a/, in 
the context of space-and-time-domain holography H20l eq. 
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(13)]. We now address the problem of optimizing the comb 
shape to maximize the efficiency. 



2. Shape optimization 

This question can be seen as a pure mathematical con- 
strained optimization. Our analysis is based on different phys- 
ical situations. The efficiency (eq. [3j is a monotonic function 
of |a_i| and ao. The largest efficiency should be obtained 
for a large |a_i | and a negative ao. The first condition makes 
sense because the retrieval is directly generated by the modu- 
lation proportional to |a_i|. Next, a negative value ao (gain) 
means that the medium is amplifying. 

This situation is realistic. One can indeed imagine that the 
atomic population is modulated between the ground state (ab- 
sorption) and the excited (population inversion). In that case, 
the absorption distribution is modulated between negative and 
positive values both bounded by au and —a m respectively 
(«m is limited by the available number of atoms). The effi- 
ciency can be larger than one. The presence of gain is interest- 
ing for classical storage by providing internal amplification. It 
has been investigated experimentally with REIC ll23ll24ll . 

The perspective of quantum storage is changing the con- 
straints. As mentioned in the introduction, the presence of 
gain reduces the fidelity of a quantum memory QS|] . Popula- 
tion inversion is then not allowed. This means that the absorp- 
tion a (A) is always positive and bounded by a^j. This is the 
situation we consider from now. 

Our analysis shows that classical and quantum storage im- 
pose different constrains. It tells us also what should be the 
shape of the comb to optimize the efficiency. Since we im- 
pose for a (A) to be always positive and bounded by au, one 
can guess that the constrains are saturated: a (A) reaches its 
boundaries. In other words, one expects that the absorption 
keeps its maximum accessible values i.e. or au- The ab- 
sorption profile should then have a square shape. The width 
of the square peaks is calculated by maximizing the efficiency 
(eq.0. 

So we define at , a 27r/T-periodic function, between —n/T 
and ir/T with a s (A) = atM between — T and T and else- 
where (see dashed lines in fig. Q]a. for example). 

rr 



The coefficients of the Fourier series are then aS 



OLM- 



and 



a M - 



sin (rr) 



r is the half-width at half- 



maximum. The proper width Tqpt that maximizes the ef- 
ficiency is 



OPT 



= — arctan 
T 



/ 2tt 



\a M L 



(4) 



At low optical thickness, Fq PT = tt/ (2T) corresponding 
to a one-half duty cycle (or finesse equals two). The width 
should then decrease for very absorbing media. The same be- 
havior has been derived by assuming that the comb is made 
of well-separated lorentzian peaks 012ill7ll . Analytic formula 
are obtained as well for the optimum width maximizing the 
efficiency (see fl2il for details). We then compare this two 
situations (see fig. [TJ. 
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FIG. 1: (color online) a. For a given maximum optical thickness 
(Qi\/L = 5 and oimL — 20 in that case), we plot the combs that 
optimize the efficiency in the two following cases: the comb is made 
of well-separated lorentzian peaks (solid black line - see 11211 for 
details) or square-shaped peaks (red dashed line), b. Optimum ef- 
ficiency for lorentzian and square shape peaks as a function of the 
maximum optical depth cxmL. The four square symbols give the 
efficiency for the shapes in fig.fjja). 

We clearly observe an improvement for the square-shaped 
peaks (red dashed line in fig. [TJa). The efficiency of the optical 
AFC is limited to 54% in the forward direction. This limita- 
tion is due to the absence of gain as discussed previously and 
is the same for the CRIB protocol [11]. This asymptotic be- 
havior is due to the trade-off between a good absorption of the 
signal and a moderate re-absorption at the retrieval. We see 
that even a limited optical thickness (typically ct]\iL ~ 10) 
gives an efficiency close to 50% for the square shape. A com- 
parable improvement is expected in the backward direction 
(Raman AFC). 



It only depends on the maximum optical thickness oimL. 
The corresponding efficiency is then given by 



?7opt — 



(XmL 



sin 2 (rg PT T) e -2aMirg PT r/7r (5) 



3. Fie/d and energy propagation 

Our analysis connects the absorption spectrum and the re- 
trieval mechanism. As previously mentioned, it may give an 
incorrect image of the protocol. This point is made explicit 
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when the energy propagation is considered in the spectral do- 
main. It is derived from the Bloch-Maxwell equations (eqs. 

®. 

d z \fl 0, oj) | 2 = -a H ® C (A) |fi (z, lu) | 2 (6) 

The optical power spectral density |f2 (z, ui) | 2 is defined by 
the Fourier transform of field in pulsed regime. The absorp- 
tion is convoluted by the homogeneous lineshape C (A) = 

7 

— — . The convolution is negligible as soon as the spec- 

7^ + A 2 

tral structures are much larger than the homogeneous width 
7. This is the case for the experiment. The equation tells us 
that the medium acts as a spectral filter. Nevertheless, it may 
be misleading to properly interpret the AFC protocol. This 
is particularly manifest for square-shaped absorption peaks. 
At large optical thickness where the highest efficiency is ex- 
pected, the peaks should be extremely narrow. The power 
spectrum is then fully transmitted (eq. [6}. In other words, 
the absorbing spectral filter has a marginal influence on the 
transmitted energy even if the protocol is very efficient in this 
regime as demonstrated in III A 21 In this situation, the en- 
ergy should mainly migrate from the signal to the first AFC 
echo (54%). The rest goes to the successive echos. This hap- 
pens without modification of the optical power spectrum. This 
property is essentially due to a frequency selective phase re- 
sponse through the propagation. As previously mentioned, the 
echo propagation obeys the causality. This assumption actu- 
ally connects the field amplitude and phase by the Kramers- 
Kronig relation. We do not explicitly calculate the phase prop- 
agation but it is embedded in the model because of the causal- 
ity. 

Our discussion is also an interesting way to underline the 
singularity of a square-shaped combs. The transmission spec- 
trum e~ a ^ L has also a square shape. The transmission 
peaks have the same width V, whatever the maximum opti- 
cal depth. This is not the case for lorentzian absorption peaks. 
In transmission, the structures become wider for a large opti- 
cal depth. They are typically broadened by a factor \[olmL- 
The transparent dispersive regions are then truncated. This 
reduces the energy transmission whose importance has been 
emphasized previously. This physical argument a posteriori 
justifies the expected improvement for square-shaped combs. 

B. Experimental square-shaped comb 

In order to support our prediction, we experimentally pre- 
pare square-shaped peaks in the absorption spectrum. We ex- 
pect a significant improvement as compared to the rudimen- 
tary preparation procedure that has been used so far (accumu- 
lation of pulse pairs see lfl6T, \vk ). This latter has the advan- 
tage of simplicity. It only needs two pulses of equal intensity. 
This produces a periodic spectrum for the optical pumping 
light. It is then not surprising that this procedure creates a 
comb like structure through the hole-burning process. The 
production of the comb has been alternatively interpreted as 
Ramsey fringes [16]. The relationship between the pumping 



light spectrum and the final comb shape is not obvious and 
is related to the atomic population dynamics. We specifically 
address this question in section III B 31 We nevertheless can 
guess intuitively that a square-shaped light spectrum with a 
2ir/T period generates a square-shaped comb through the op- 
tical pumping process. We now prepare this specific optical 
spectrum and use it for the preparation procedure. 

1. Square-shaped optical spectrum with a 2ir/T period 

A 27r/T-periodic spectrum for the pumping light \£ (u) | 2 
is produced by a pulse train in the time-domain. The tem- 
poral separation T imposes the period and the pulse duration 
gives the total bandwidth. The fine structure of the spectrum 
depends on the relative pulse amplitudes. 

Experimentally, the temporal shaping of the preparation 
pulses is realized by an acousto-optic modulator (AOM) con- 
nected to an Arbitrary Waveform Generator (AWG 520 Tek- 
tronik). We here use the same optical setup that have been 
previously described (17]]. This allows a complete control in 
amplitude and phase of the preparation pulse sequence £ (t). 
We'd like to obtain a a field spectrum £ (u>) which has a con- 
stant value between — T and T and elsewhere in the interval 
[— tt/T; tt/T]. In the time-domain, it is a pulse train with rel- 
ative amplitudes P (k) — sinc(fcrT/2) (k is integer, P (k) 
can be negative). 

By controlling the AWG RF amplitude and phase sent to 
the AOM, we create an arbitrary sequence of pulses. The non- 
linear RF response of the AOM should be taken into account 
in order to faithfully control the field spectrum. The RF to op- 
tical intensity response is recorded and is used as a correction 
factor to obtain the P (k) amplitudes of each pulse (for neg- 
ative values of P (k), the RF phase is flipped). We truncate 
the series to —30 < k < 30 corresponding to 61 pulses. The 
pulse duration is 300ns adapted to the AOM bandwidth and 
the temporal separation is 1.5/is (1/T = 666 kHz) ifnll . We 
experimentally verify the field spectrum by heterodyning the 
pulse sequence with a shifted reference from the same laser. 

We observe in fig. |2]a good overall control of the field spec- 
trum £ (u>). By applying 61 pulses, we obtain a square-shaped 
spectrum with a 2tt/T period. The peak width is changed by 
modifying the pulse sequence (see inset). 

The spectrum is asymmetric about the carrier frequency. 
This is explained by an asymmetry in the AOM frequency re- 
sponse. It strongly depends on the alignment of the device 
itself and the optical fiber placed behind. More importantly, 
we note an imperfect contrast in the modulated spectrum. This 
is especially true for its central part. It can be due to a rela- 
tive distortion of the pulses. This point is not surprising be- 
cause of the previously mentioned non-linear response of the 
AOM. The potential effect of the pumping dynamics will be 
discussed later (see section III B 3\ . We prepare four differ- 
ent time sequences. The first is a pulse pair (PP-sequence) 
and is taken as reference. The three others S1/2, S1/3 and 
S ^-sequences have a square-shaped spectrum with different 
widths respectively tt/ (2T), tt/ (3T) and tt/ (5T). We now 
experimentally evaluate their influence on the AFC efficiency. 



5 




Frequency (MHz 



FIG. 2: (color online) Heterodyne beating of the preparation pulse 
sequence for T = it/ (2T) with l/T = 666 kHz (red dashed line). 
Since the reference is monochromatic, it directly gives the field spec- 
trum \£ (w) |. We compare it with the spectrum of a pulse pair 
P(0) = P(l) = 1 (solid black line). Inset: Central part of the 
spectrum [—l/T; l/T]. We add the pulse sequence corresponding 
to r = 7r/ (5T) that is also used in the experiment (magenta dashed 
dotted line). 



2. Measured AFC efficiencies 



We implement the protocol in 0.5% doped-Tm 3+ :YAG 
crystal [17]. A 210G magnetic field is applied along the [001] 
crystalline axis and splits the ground and excited levels into a 
nuclear spin doublet by A s = 6MHZ and A e = 1.3 MHz re- 
spectively. The comb structure is obtained by optical pumping 
between these long-lived spin levels. The laser polarization is 
also applied along the [001] axis. The crystal is immersed in 
liquid helium at 2.3K. This ensures a narrow homogeneous 
linewidth (7 ~ 2ir x 5kHz) and then an accurate spectral tai- 
loring. 

The four previously described pulse sequences (PP, S1/2, 
Si/ 3 and Sx/5) represents an elementary pattern whose total 
duration is 100/is. This pattern is repeated 5000 times to ac- 
cumulate the population in the shelving state with a low pump- 
ing power (weak area pulses). The resulting comb depends on 
the pulse sequence and the pumping power. It is probed after 
50ms waiting time for a complete decay of the excited state 
population. 

The independent weak probe beam is controlled by a sec- 
ond AOM for frequency scanning and temporal shaping. First, 
it is used to measure the comb shape by sweeping the AOM 
carrier frequency. We obtain the absorption spectrum (see 
fig [3] a. and c). Next, a 450ns signal pulse is sent into 
the medium. It is longer than the preparation pulses to en- 
sure a good spectral overlap of the signal bandwidth and the 
comb. This probe pulse is absorbed and retrieved at a later 
time T = l.b^ts (see fig[3]b and d). 

We see a net gain between the PP and the Si/2-sequence 
since they give an efficiency of 10.1% and 17.3% respectively. 
This is expected since the two sequences generate different 
comb shape (see fig[3]a and c). The sharp edges of absorption 
spectrum approaching a square shaped structure for the S1/2- 
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FIG. 3: (color online) a. and c. Central part of the absorption spec- 
trum after the PP and S1/2 sequences respectively. We fit the spectra 
by a simple optical pumping model (dashed blue line, see section 
III B 3l for details). The insets show the same spectrum on a larger 
bandwidth, b) Optical AFC echo of a weak 450ns signal pulse af- 
ter the PP-sequence (corresponding to a), d. Optical AFC after the 
Si/2-sequence (corresponding to d). The pumping power is adjusted 
independently for both sequences to optimize the retrieval efficiency. 



sequence (see fig [3] c) explains a gain in efficiency. There 
are also noticeable differences between the field spectrum and 
atomic comb (for example the inset of fig. [2] and the fig. 
[3J: can be compared). First, the absorption spectrum is vis- 
ibly smoothed out by the optical pumping dynamics. More- 
over, we cannot obtain a fully contrasted comb. More pre- 
cisely, the transmission should be 1 at the maximum and few 
percent at the minimum limited by the sample optical depth 
(cxmL ~ 4 — 5 in our case). The maxima and minima of the 
absorption spectrum strongly depend on the pumping power. 
The latter is chosen to optimize the retrieval efficiency for the 
PP and S1/2 sequences. In order to have a good understanding 
of the preparation population dynamics, we vary the pumping 
power and record the corresponding efficiency (see fig. |4|. 
We decide to use the mean transmission value (central part 
of the absorption spectrum) as abscissa to plot the different 
efficiencies. It as the advantage to be easily accessible exper- 
imentally. It is also obtained from the model that we develop 
in section HI B 31 

We expect different results for different square widths as 
previously discussed in section QTA~2] We indeed see that the 
square width should be properly chosen to obtain the high- 
est efficiency (eq. @). We then also use the sequences S1/3 
and S1/5 to observe this dependency. We recognize the fea- 
tures derived from our analysis of the optimum comb shape 
(section Til A 21 ) i) a square like shape of the comb improves 
significantly the retrieval, ii) the efficiency strongly depends 
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Average comb transmission 

FIG. 4: (color online) Efficiency comparison between the PP, S1/2, 
1/3 and S1/5 sequences. The lines connecting the measurements are 
used to guide the eye. The four insets show an absorption spectrum 
for each sequence. The corresponding measurements are indicated 
by arrows. The two circles represent the measured efficiencies and 
spectra that have been analyzed in fig[3] 



on the square width of the optical spectrum. 

We should be able to predict the efficiency from the atomic 
comb shape (eq. [3). This procedure demands a fine knowl- 
edge of the absorption spectrum i.e. the average value a>o and 
the modulation a_i. We measure both coefficients by record- 
ing the transmission spectrum. The results of eq. [3] can then 
be compared to the independently measured efficiencies. This 
approach is unfortunately unsuccessful to fit the data set pre- 
sented in fig. H] To explain this discordance, we see two rea- 
sons. The first one is technical. It is not easy to accurately 
measure a large optical thickness. In that case, a weak trans- 
mission can be biased by experimental artefacts [17]. The 
second one is fundamental. We implicitly assume in section 
III A 1 I that the atomic comb is uniform along the crystal length. 
This allows a complete analytical calculation. In the contrary 
case, «o an d ot-i (in eqs. |2]i would be z-depend. We can ques- 
tion this assumption. The comb is produced by a specific pulse 
preparation sequence. During the propagation, this sequence 
can be distorded by the atomic comb under construction. As a 
consequence, the resulting comb would not be uniform along 
the crystal. In order to analyze the impact of these two effects 
and have a better physical understanding of the AFC limita- 
tion, we develop a simple pumping model to discuss the pre- 
vious data. 



3. Discussion - optical pumping dynamics 

The optical pumping dynamics is a key element to have at 
least a qualitative understanding of the current limitations. We 
here briefly develop a model to account for the frequency se- 



lective hole-burning preparation. 

The rate equation limit of the Block-Maxwell equations 
have been successful to predict the hole shape in high reso- 
lution spectroscopy with monochromatic pumping [21]. This 
approach is perfectly valid in our case since each sequence 
pattern satisfy the weak area approximation. In that case, 
each sequence (duration T p = lOOfis) induces an incremen- 
tal small population change w — > w + 5w that is given by 

5w = —C (A) (g) |£ (A) \ 2 w. The optical spectrum |£ (A) | 2 

(£ is the Fourier transform of the pulse sequence) is convo- 



luted by the homogeneous lineshape £ (A) 
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The 



population difference between the excited and the ground state 
is w — n\ — n 2 . 

Since we repeat each sequence to accumulate the popula- 
tion in a shelving state, we model the preparation by a contin- 

uous frequency selective pumping rate R (A) 



wT„ 



R(A) 



2-kT, 



■C (A) ® \£ (A) 



(7) 



The hole-burning process is due to at least one shelving state. 
In our case, even if most of the population ends up in the long- 
lived Zeeman state, the pumping dynamics also involve the 
3 F4 bottle-neck state (lifetime 10ms). To simplify the prob- 
lem, we model our system as an equivalent three level system 
(see fig. 01. 
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FIG. 5: Left: Thulium level structure. Right: Equivalent three-level 
system for optical pumping model. 

The rate equations for the populations ni,n 2 , and n 3 (see 
fig. |5]l are then 

d t ni (A) = \R (A) (n 2 - n x ) + ^n 2 - {n x - n 3 ) 

d t n 2 (A) = ±R(A) (m - n 2 ) - ±-n 2 



d t n 3 (A) 



2 

l-r „ 
Ti 



-n 2 



h ("-3 - ni) 



Tz 



(8) 

At the end of the pumping sequence, the populations reach 
a steady state nf,n 2 , and n§. Before sending the signal, 
we wait for a complete decay from the excited state. The 



populations should now be n x 



a P 

rn 2 ,n 2 



0, and 



13 = 7J3 + (1 — r) n 2 . The absorption spectrum on the 1-2 
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transition is given by 



a (A) 



where e = 



a M nf (A) 



1 



OiM- 



{l + r)R(A)T 1 + 2 
R{A)T 1 /e + 4 



(9) 



. The parameter e should be 

[l-r)T z /Ti + 3 
much smaller than one. The population then accumulates is 3 
because (1 - r) /T x » l/T z . 

This formula connects the optical spectrum and the absorp- 
tion profile. To obtain an analytical solution for the AFC ef- 
ficiency, we have to assume that the optical spectrum is uni- 
form over the entire length of the medium. In that case only, 
the atomic comb is uniform as well. That is an underlying as- 
sumption of the efficiency formula (eq. |3)- This supposition 
will be discussed later on. 

Such a model allows us first to predict the absorption spec- 
trum for any predefined pumping sequence. We here focus on 
the PP and Si/2-sequences which yield the highest efficien- 
cies. The PP-sequence produces a cosine-modulated spec- 
trum. In our case, since 1/7 >> T, we obtain for the pump- 
ing rate at the center of the spectrum (eq. |7): R (A) oc 
1 + exp (7T) cos (AT). By using the formula [9] we fit the 
absorption spectra produced by the PP-sequence with differ- 
ent pumping powers. We see a typical comparison between 
our model (dashed blue line in fig[3]a) and the experiment. 

The Si/2-sequence has been designed on purpose to gener- 
ate a square-shaped optical spectrum. The convolution by a 
lorentzian smooths out the edges and the pumping rate now 
,v-.__.._/'27r(p+l/2)-Ar i We 



reads as R (A) oc 2 arctan 



7 T 



also fit the experimental absorption spectrum in fig[3]c). 

The agreement is satisfying for both sequences. Our model 
is able to properly predict the observed transmission spectrum. 

We now interpret the connection between the heterodyne 
beat and the absorption spectrum. We have seen that the 
Si/2-sequence induces a residual background reducing the 
contrast of the modulated spectrum (see fig. |2). Our model 
does not include this component but properly reproduces the 
observed absorption spectrum. We conclude that this back- 
ground does not have a major impact. It seems that the reduc- 
tion of the comb contrast through the optical pumping dynam- 
ics is mainly due to the homogeneous linewidth (convolution 
effect). This limitation is fundamental and not technical. 

The comparison with the observed absorption spectrum val- 
idates our approach of modelling the pumping process. Com- 
bining eqs. [9] and [3] we are able to predict the observed ef- 
ficiency (in fig. |4|i for the PP and Si/2-sequences. As we do 
not independently measure the pumping power which gives 
the amplitude of £, we also use the mean transmission value 
(central part of the absorption spectrum). The average trans- 
mission is indeed predicted by our model (eq. [9]) and is used 
in abscissa. We then plot the efficiency for different pumping 
power (see fig. [6]). The variable pumping power parametrizes 
the curve. The only remaining adjusting parameter is the ini- 
tial optical depth oimL. 

The initial absorption is typically olj^L ~ 4 — 5. It is not 
completely easy to accurately measure a large optical depth. 




0.2 0.4 0.6 0.8 1 

Average comb transmission 

FIG. 6: (color online) Comparison between the measured (same as 
in fig. [4]l and the predicted efficiencies for the PP (a) and S1/2- 
sequences (b). For each sequence, we plot the prediction of our 
pumping model with cimL — 3 olmL = 6 for the PP (a - solid 
black lines) and Si/ 2 -sequences (b - solid red lines). 



The efficiency strongly depends on a^jL. To illustrate this 
dependency and consider a larger range, we plot the expected 
efficiencies for ojhL = 3 and a^jL = 6 (solid black lines 
and red dashed lines in fig. [6] a and b for the PP and S1/2- 
sequences respectively). 

The general qualitative behavior is properly reproduced. 
Without pumping light, the mean transmission is e~ au and 
the efficiency zero (the comb is absent). At very high pump- 
ing power, the medium is completely transparent and there is 
no AFC echo neither. The power broadening effect is indeed 
smoothing out the periodic modulation of the pumping spec- 
trum and uniformly bleaching the sample. Quantitatively we 
observe a major discrepancy between the curves. We cannot 
reproduce the position of the maximum even if we consider 
optical depths (qj/ L — 3 and 6) which are largely surround- 
ing the experimental value. This may sound surprising since 
we fit accurately the measured transmission spectrum (see fig 
[3] a and c). To explain this discordance, we go back to the 
conclusion of section llIB2l If the atomic comb was uniform 
over the entire length of the medium, we should be able to 
fit both the measured transmission spectrum and experimen- 
tal efficiencies. Our optical pumping model predicts indepen- 
dently both and should allow us to get around the experimental 
artefacts biasing the transmission measurements. We then ex- 
pect that the atomic comb is actually non-uniform. This would 
certainly strongly modify the propagation of the AFC echo (in 
eqs. |2j ao and a_i would be z-depend) and change the final 
efficiency. But it sounds possible to fit properly the finally 
observed transmission spectrum. It indeed results from an av- 
erage absorption spectrum over the propagation but could be 
fitted by an equivalent uniform profile. In that sense, the AFC 
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efficiency would be more sensitive to the uniformity of the 
comb than the transmission spectrum. 

The propagation of the pumping spectrum through the sam- 
ple produces certainly a non-uniform atomic comb. To ac- 
count for this propagation effect, one can include the back- 
action of the absorption spectrum (propagation of the pump- 
ing light). This leads to two coupled equations that can be 
solve numerically. This work is currently under investigation 
in our group. 



III. INFLUENCE OF THE MAGNETIC FIELD 

The previous analysis is based on an optical pumping 
model. We simplify the level structure (see fig. [3j to obtain an 
equivalent structure and stay close to the physical intuition. If 
we neglect the influence of the 3 F4 bottle-neck state, such a 
simplification seems to be appropriate as soon as the ground 
and excited state Zeeman splittings (A g and A e ) are much 
larger than the comb bandwidth. The shelving state is then 
out of resonance. Spectral hole-burning (SHB) is an alterna- 
tive way to consider the pumping dynamics. SHB has been 
widely used as a fine spectroscopic tool with monochromatic 
laser especially for Zeeman interaction and superhyperfine in- 
teraction 12511 . Under magnetic field, the spectrum is usually 
enriched by side-structures as holes or anti-holes. Following 
this analogy, the preparation of a comb by optical pumping 
should also lead to side combs and anti-combs. If this struc- 
tures are completely out off resonance, one could indeed treat 
the preparation as a population accumulation in a single far- 
off-resonance shelving state. For a weak magnetic field or 
when a large bandwidth is demanded, the side structures over- 
lap the main comb. They can be caused by the Zeeman (hy- 
perfine) or the superhyperfine effect. We specifically study 
these situations here. 



A. Zeeman effect 

The magnetic field is applied along the [001] crystalline 
axis. In that case, all the excited sites are equivalent. The 
Zeeman splittings have been measured previously: 5 g = 28 
kHz/G and 5 e — 6.0 kHz/G respectively for the ground and 
the excited state [263. The experimental 210G magnetic field 
produces shifts of A g = 5.9 MHz and A e = 1.3 MHz. 
The AFC bandwidth is typically 4MHz limited by the 300 ns 
preparation pulse duration (see inset of fig. [3] c. for exam- 
ple). Even if the ground state splitting is larger than the band- 
width, this not the case for the excited state. This modifies the 
pumping dynamics. In SHB spectroscopy a monochromatic 
laser creates a hole at the central frequency. In the present 
situation with four levels, we observe also two side holes po- 
sitioned at ±A e , four anti-holes at ±A S and ± (A ff — A e ). 
Our preparation procedure is then creating a comb centered at 
the central frequency, two side combs at ±A e and four anti- 
combs at ±A g and ± (A s — A e ). Since A g >> A e , the 
influence of the anti-combs should be negligible for a 210G 
magnetic field. Nevertheless we cannot neglect the two side 



combs at ±A e = ±1.3 MHz as it is smaller than the band- 
width. To minimize the interplay between these structures, 
we have intentionally chosen the magnetic field in order to 
have A e ~ 2/T. This a priori assumption deserves further 
study. To question it, we use the PP-sequence. We then sim- 
ply measure the optical AFC efficiency as a function of the 
applied magnetic field (see fig. [7). 
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Magnetic field (Gauss) 



FIG. 7: AFC efficiency as a function of the magnetic field (the solid 
line is used to guide the eye). The error bar are smaller that the 
marker size. The 210G value that has been used for the previous 
experiments is indicated by an arrow (see text for details). 

We observe a strong influence of the magnetic field on the 
efficiency. We analyze this figure as follows. 

At low magnetic field (<50G), the efficiency is zero. This is 
directly due to the population lifetime in the Zeeman state. We 
cannot observe long-lived holes by SHB spectroscopy. The 
preparation is then inefficient. We observe an echo only if the 
population lifetime is significant (> 100G in our case). 

The range 100- HOG is very singular since A g and A e are 
both smaller than the bandwidth (~ 4MHz). The situation 
is complex because the main comb, the side combs and anti- 
combs overlap each other. But nevertheless it leads to a high 
efficiency. If the side combs are shifted by the main comb 
spacing (A e = p/T, p is integer), and the anti-comb by half 
of it (A 9 = (pi + 1/2) /T, pi is integer), the resulting pro- 
file is very contrasted. In other words, the preparation spec- 
trum is optically pumping and de-pumping at the same time 
within the four-level system. The two mechanisms act ac- 
cordingly to increase the comb contrast. We cannot satisfy 
all the conditions experimentally when A g and A e are kept 
smaller than the bandwidth. A 130G magnetic field repre- 
sents the best trade-off between the conditions: A e = 1.2/T 
and A g = 5.5/T. We indeed observe a local maximum of the 
efficiency for this value (see fig. |7). It supports our analysis. 

In the range 140-430G, only A e is smaller than the band- 
width. An overlap of the side combs with the main one is 
possible. The matching is realized for 220G (A e = 2/T) and 
330G (A e = 3/T). We also observe a local maximum of the 
efficiency at 220G close the value chosen for the previous ex- 
periments. For higher values of the magnetic field, the overlap 
is decreasing and the efficiency constant. 

Our study is particularly relevant for the optical AFC. We 
clearly see that high efficiencies are obtained when the four- 
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levels are all contained in the interaction bandwidth. Such 
a configuration is less suitable for the Raman AFC. First, a 
spin doublet in the ground state is not sufficient to provide one 
shelving state plus two spin states (long coherence lifetime). 
Next, the Raman transfer demands to address independently 
the different transition. A overlap may not be appropriate es- 
pecially for the ground state where the coherences are stored. 
A partial overlap in the excited state may be tolerated for the 
Raman AFC but this possibility has not been investigated so 
far. 



B. Superhyperfine effect 

The superhyperfine effect corresponds to the next order of 
magnetic field induced interaction. It is produced by the sur- 
rounding nuclear spins on the dopant rare-earth l27ll . It is usu- 
ally studied with EPR or NMR spectroscopy. An optical de- 
tection makes the observation easier for luminescent centers 
ll25l p.72]. SHB spectroscopy can also be used to study this 
ion-ligand interaction [28]. The superhyperfine effect is then 
indicated by the presence of adjacent anti-holes on the broader 
hyperfine or Zeeman structure. We perform SHB measure- 
ments for a varying magnetic field along the [001] crystalline 
axis (see fig. ©. 
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FIG. 8: (color online) Spectral hole-burning spectroscopy of the Zee- 
man excited state splitting (blue markers) and the superhyperfine in- 
teraction (red markers). We plot the side holes position as a function 
of the magnetic field. Insets: Examples of the hole-burning spectra 
(for 1 10, 210 and 580G from left to right). We confirm previous mea- 
surements for the excited state Zeeman shift S e — 5.91 ± 1% kHz/G 
12611 . For large magnetic fields (>210G), we see adjacent side holes 
that we attribute to the superhyperfine interaction with Al nuclei. The 
splitting is 5s = 1.05 ± 1% kHz/G. 

We initially observe side holes at ±A e due to the Zee- 
man interaction (see fig. [8}. The previously mentioned anti- 
holes at ±A g and ± (A 9 — A e ) are present but rapidly out 
of the spectrum range and small because of the large optical 
depth. The strongest surrounding nuclear moments are the 
aluminium ones. In that case, the population distribution is 



usually governed by Am/ = ±1 spin flips [29]. Our result 
6s = 1-05 ± 1% kHz/G is consistent with the Al gyromag- 
netic ratio 1.11 kHz/G. One can be surprised by the intensity 
of the superhyperfine effect (hole depth) which is comparable 
to the Zeeman one. Nevertheless the Al ions occupy differ- 
ent sites in the crystal leading to different coupling strengths. 
Since we do not observe any quadrupole structure, we cannot 
distinguish them [30]. Such an analysis demands an accurate 
spectroscopic study. 

The superhyperfine interaction modifies the AFC prepara- 
tion dynamics and generating side structures as the Zeeman 
effect does at lower magnetic field. At moderate magnetic 
field, the shift is smaller than the AFC peak and then broaden 
it. This reduces the contrast of the comb and then the ef- 
ficiency. In our case, the peak width is typically 200kHz 
corresponding to a 180G superhyperfine shift. For a larger 
magnetic field, side structures should clearly appear. Even if 
because of the high number of side structures (Zeeman and 
superhyperfine), the matching between them is relatively un- 
likely for a given magnetic field. 

The effect should be stronger for Kramers ions in high spin 
concentration matrices. Its influence on the AFC storage has 
already been observed 11611 . A systemic material study should 
be pursued especially because a high bandwidth and then a 
large magnetic field is desired for quantum repeaters applica- 
tion. 



IV. CONCLUSION 

We carefully study the efficiency of the optical AFC pro- 
tocol. We show that the preparation procedure can be sig- 
nificantly improved. The most efficient atomic comb have 
a square shape whose width should properly adjusted for a 
given optical depth. This requires an accurate control of 
preparation pulse sequence in amplitude and phase. We im- 
plement this technique in a Tm 3+ :YAG crystal and increase 
the efficiency from 10.1% to 17.3%. The experimental results 
are consistent with our expectations. We develop a simple 
rate equation model to describe the optical pumping dynam- 
ics. Even if it properly reproduces the comb spectrum, we are 
unable to predict ab initio the observed efficiencies. This dis- 
crepancy tells us that the propagation of the preparation pulses 
inside the medium and the backaction of the prepared atomic 
comb should be taken into account. This route is under inves- 
tigation in our group. 

We also study the effect of the magnetic field. This is rel- 
evant in the perspective of high bandwidth application. The 
Zeeman and the superhyperfine interactions have a strong in- 
fluence on the optical AFC protocol. They can be beneficial 
or detrimental. We clarify these conditions. 

Our work should be placed in the context of quantum stor- 
age where a Raman conversion has to be implemented (Raman 
AFC). Our analysis mainly focuses on the preparation stage of 
optical AFC but can be extended to understand the efficiency 
of the complete protocol. 
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